Vectorial Phase Retrieval for Linear Characterization of Attosecond Pulses 
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The waveforms of attosecond pulses produced by high-harmonic generation carry information 
on the electronic structure and dynamics in atomic and molecular systems. Current methods for 
the temporal characterization of such pulses have limited sensitivity and impose significant ex- 
perimental complexity. We propose a new linear and all-optical method inspired by widely-used 
multi-dimensional phase retrieval algorithms. Our new scheme is based on the spectral measure- 
ment of two attosecond sources and their interference. As an example, we focus on the case of 
spectral polarization measurements of attosecond pulses, relying on their most fundamental prop- 
erty - being well confined in time. We demonstrate this method numerically reconstructing the 
temporal profiles of attosecond pulses generated from aligned CO2 molecules. 

PACS numbers: 42.30.Rx,78.47.D- 
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Optical and XUV pulses with durations significantly 
below one picosecond cannot be directly characterized in 
the time domain, since there are no suitable photodetec- 
tors. Therefore they are usually characterized in the fre- 
quency domain by measuring both the spectral amplitude 
and spectral phase of the pulse. The former may be de- 
termined straightforwardly by means of a spectrometer. 
The latter requires either a fast modulator, a detector 
or a reference pulse with which the unknown pulse can 
be interfered. A modulator or gate of sufficiently rapid 
response is typically synthesized by means of nonlinear 
optical processes In the femtosecond regime, most 
measurement schemes are based on nonlinear light mat- 
ter interactions. By exploiting media having a nonlinear 
response, various successful schemes have been developed 
for complete characterization of femtosecond pulses, most 
notably FROG and SPIDER 0. 

The ability to produce attosecond pulses has set new 
benchmarks in time-resolved measurements. Such pulses 
enable one to probe electron dynamics on the atomic time 
scale. Recently it has been demonstrated that the at- 
tosecond production process carries information about 
both electron dynamics [1] and molecular structure 
However, full access to the information contained in the 
pulse requires its complete characterization. A direct im- 
plementation of the pulse characterization schemes de- 
veloped for the femtosecond regime is challenging, due 
to the low signal levels and the absence of appropriate 
nonlinear media in the XUV. 

Several characterization schemes have been devel- 
oped for the attosecond XUV domain How- 
ever, time resolved measurement of attosecond pulses re- 
mains a major challenge. For example, certain impor- 
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tant attosecond-scale processes such as plasma mirrors 
(1^ and field-enhancement HHG [3l have not yet been 
fully characterized. 

In this letter we propose a new approach for measuring 
attosecond pulses, which - in contrast to other methods, 
relies only on linear spectral measurements of the radia- 
tion generated by several sources. The key new feature of 
the method that enables this approach is the explicit uti- 
lization of a temporal support constraint in the retrieval 
algorithm. The fact that the pulse is limited in duration 
allows a unique solution to be extracted. To illustrate the 
approach, we consider the two polarization components 
of the attosecond pulse as two independent sources. We 
show that the spectral measurement of the two polar- 
izations together with their relative phases, obtained by 
means of spectral interference, is sufficient to retrieve the 
spectral phase of each component individually. We nu- 
merically test the method on attosecond pulses calculated 
for aligned CO2 molecules , in which the polarization 
varies non-trivially with the frequency. We establish ex- 
cellent reconstructions even for pulses that do not exactly 
satisfy the finite duration assumption, but rather have a 
sufficiently rapid decay. 

Phase retrieval problems are common in many 
branches of physics, including astronomy 15| , NMR UJ, 
ultra-fast optics [2|, crystallography [13 and lens- less 
imaging [18|-|2(J|. In each of these cases, the measured 
signal spectrum, together with some assumption on the 
signal itself, allows one to retrieve the phase, without di- 
rectly measuring it. In our case, the assumed property is 
that the pulse has a limited, but not necessarily known, 
"time window" in which the intensity is non-zero, usually 
referred to as compact support. Generic compact support 
phase retrieval problems are known to have many solu- 
tions in the one-dimensional (ID) case, but usually a 
unique solution for higher dimensions [21|. As our mea- 
surements are of a one dimensional Fourier transform, the 
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compact support assumption is not enough to reconstruct 
the spectral phase. This is solved by measuring the rel- 
ative phase between two (or more) spectra. Such a mea- 
surement enables a vectorial ID phase retrieval problem. 
If the two components of the vector are non-degenerate 
(in a sense to be defined later on), then the problem has 
a unique solution for the spectral phase, up to a phase 
linear in frequency (i.e. a arbitrary delay which is not 
physically significant). It is useful to consider how the 
temporal support constraint enables a unique inversion 
of the data. We first develop a formal approach that 
illustrates this, and then describe a simpler inversion al- 
gorithm. 

Consider the electric field spectrum E{uj) of a pulse 
with a finite duration T, sampled p8| at coj = 27r/T 
for j ~ 1,...,N. The spectrum can be written as 
Eiu}) = Et E{t)e~''^* = Y.t E{t)z\ where z = e"*'^. Ac- 
cording to the fundamental theorem of algebra, we can 
write E{z) ~ E{t)z* = Yij ~ ^j) where Zj are the 

roots of the polynomial E{z). A linear measurement 
of the spectrum measures |i?(a;)p, which at the sampled 
points l^l = 1 can be shown to equal: 

1^(^)1' = 7^ n - ^^^(^ - ^-') (1) 

where is the complex conjugate of Zj. Without any 
prior knowledge about the pulse, the spectral phase is, 
by definition, unrecoverable. This is because the polyno- 
mial representing the spectrum is under-sampled: E{uj) 
corresponds to a polynomial of degree A^ in z, but |i?(a;) p 
corresponds to a polynomial of degree 2A^ in z. There- 
fore, the A^ samples of [^^(z)^ are insufficient to unam- 
biguously determine E{t). If, however, a compact sup- 
port constraint is assumed, meaning E(t) = for, say 
the A^/2 samples in the range T/2 < t < T, the degree 
of |£'(z)p is at most A^. Therefore, it is well sampled 
and the A^ roots of the polynomial (zj and zj"^) are 
uniquely determined. The compact support constraint, 
nevertheless, does not remove all the ambiguities. In or- 
der to retrieve the pulse E{z), only one root from each 
pair of roots, zj and 'Zj~^, should be chosen. As there 
are 2^/^ possibilities to choose one root from each pairs 
of |i?(2;)^|, there still exist 2^/^ diff'erent pulses, all hav- 
ing both the measured spectrum and the correct compact 
support. This ambiguity in the ID phase retrieval prob- 
lem is well known [2l|. Hence, additional information is 
required to unequivocally determine the "correct" choice. 

The phase ambiguity can be overcome by using two 
(or more) spectra which we refer to as "components" , 
and the relative phase between them. We will de- 
note such measurements as "vectorial" . As an exam- 
ple of a vectorial measurement, appropriate for the at- 
tosecond domain, we use polarization: when attosecond 
pulses are generated from an anisotropic media, such 
as aligned molecules by means of high-harmonic genera- 
tion, non-trivial frequency-dependent polarization is ex- 
pected [12]. For the two component spectra, |i?:c(aj)p 



and \Ey{uj)\'^ with the same time domain compact sup- 
port, we can find the A^ roots {zx}j, {zy}j correspond- 
ing to Eq. dD) and their complex conjugates. Without 
the relative phase, these are two independent 11? phase 
retrieval problems, each having many solutions. How- 
ever, a complete polarization measurement, consisting of 
both the spectra of two orthogonal polarizations as well 
as the interference spectrum between them, provides suf- 
ficient information to eliminate the ambiguities. Not- 
ing that Ex^y{Lu) = \Ex,y{uj)\e'^''-^^'^\ the vectorial mea- 
surements allows us to establish \Ex{uj)\'^, \Ey{uj)\'^ and 
Ex{u)Ey{oj) = \E.x\\Ey\e''-'t'^-'^y\ These quantities are 
represented by the factored polynomials: 

mz)\' = ^% Uiz- z^){z - If"') (2) 

^ > 3 

\Ey{z)? = n - - (3) 

E.Az)Ejz) ^ ^HiL (^ „ ^^-)(^ _ ■^-') (4) 
j 

The essential point is that from these equations the 
"which root" ambiguities of both components can be re- 
solved: the correct roots for Ex{z) are those that are 
common to |i?2:(z)P and Ex{z)Ey{z), and similarly for 
Ey{z). In Figlll we show an example of the method for 
a simple pulse. The time domain signals are shown in 
the inset. In the main figure, the roots generated by 
Eqs. (|2l3p are shown as dots and crosses in the complex 
plane (the black line is the unit circle). As expected, the 
roots come in pairs: for example, the roots marked by A 
and B are related by 2 a = 'zb~^ ■ The one dimensional 
phase retrieval problem requires one to choose the cor- 
rect root from each pair. In the vectorial case, this can 
be done using the roots generated by Eq.(|3]): these roots 
are marked by green squares. As seen in the figure, each 
square coincides with only one of the roots, thus identi- 
fying the correct root of each pair. In the above example, 
A, rather than B, is the correct root, as it coincides with 
a root of Eq. (g]). 

Clearly, if |i?2,(z)p and |i?j,(z)^| have m common ze- 
ros, there is still an ambiguity among 2™ different pulse 
shapes. We will call such cases degenerate. In the absence 
of noise, the set of pulses which share common zeros is 
very small, unless they have been manipulated in a com- 
mon manner to introduce them. Nevertheless, common 
zeros do arise when both pulses are passed through the 
same linear phase-only filter. The pulses then share a 
component of common spectral phase, leading to degen- 
eracy. In this case the pulse fields cannot be retrieved 
by the above method. Further, most experimental data 
will have some noise S{z). This will tend to move the 
zeroes of the polynomials from Zi to Zi — IS.Zi . A simple 
calculation shows that 

E{z) + S{z) = E{z) + Az, Y[{z Zj) + 0(Az2) 
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FIG. 1: An example for our reconstruction procedure for two 
artificially generated pulses with the same compact support. 
The blue dots, red crosses and green squares are the roots, 
in the complex plane, of |_E2,(2)|'^,|i?j,(2:)|'^ and Ex{z)Ey{z) 
respectively. A time domain plot of the two polarizations is 
shown in the inset. 



Assuming that all the zeros are shifted by about the same 
distance , the position uncertainty radius scales with N 
as Az ~ In contrast, the distance between the zeros 

scales with as 1/\/N. Hence, increasing the spectral 
bandwidth (and therefore increasing N without changing 
the compact support) allows the zeros to be completely 
distinguished. 

Another consideration is that most physical pulses do 
not have exact compact support, but rather a sharp decay 
outside some temporal region. Our method is applicable 
even in such cases, as we now show by means of an alter- 
native solution method. Ex{t) and Ey(t) can be written 
as: 

E^it) = ^i?,(L.)e'"* = J2 lEA^W^'Xiiu) (5) 
Ey{t) = |£;.y(a;)|e'"*e'('^«(")-'^-("»X(w) (6) 

where X{u;) = e^'t'^i'^) represents the spectral 
phase. Applying the compact support assumption gives 
E^x{t)\t-i N = Ey[t)\^^^ jv = 0. Using this condition 
in Eq. (|5l6p gives a set of N linear homogenous equations 
for the A'^ unknowns X{ijj). Exploiting the arbitrariness 
of the absolute phase, one can set X{u:i) = 1 and obtain 
an over-determined set of inhomogeneous equations. By 
solving this set, one can solve the phase problem provided 
X{ui) is a phase-only function, that is, |^(a;)| = 1 for all 
Lo. This has two important consequences: (1) Since the 
problem boils down to solving linear equations, the sen- 
sitivity to noise is linear in the noise amplitude; (2) In 
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FIG. 2: (a) Two time domain polarization components of the 
attosecond pulse generated from CO2 molecules at an align- 
ment of 58° (between the molecule symmetry axis and the 
IR polarization) and IR laser intensity of 0.07aM at 800nm. 
(b) Relative change in the pulse energy ( defined in Eq[7ll for 
different assumed compact supports. 



cases where the solution to Eq. ()5|6|) is far from yielding 
a unimodular complex number, the compact support as- 
sumption is either wrong or there is a degeneracy. As we 
later show, this gives us the means to search for the cor- 
rect compact support without assuming it beforehand. 

Our proposed algorithm seeks a domain of compact 
support choosing that which is most consistent with the 
above conditions as our constraint. For each assumed do- 
main of support, r, we resample the spectral information 
at discrete frequencies with spacing VI = ^-k/T. We then 
use £^2-,y(i)lt=i,...,.N. = in Eqs. (O and ^ to find X{uj). 
For each T we calculate how far X{u:) are from being 
unimodular complex, by calculating the relative change 
in the pulse's energy when using X{uj) (which might not 
be unimodular) as the spectral phase: 



Err{T) 



EJE{uj)m-X{u;U^ 



(7) 



Err{T) is the standard metric used in the Gerschberg- 
Saxton algorithm [15[ for compact support phase re- 
trieval problems. A plot of Err{T) is given in Figl5] 
We then choose the value of T that minimizes Err{T) 
as the compact support domain and the argument of the 
corresponding Xiuj) as the pulse's phase. We demon- 
strate our method by reconstructing simulated attosec- 
ond pulses generated from aligned CO2 molecules (for 
HHG simulation details see 141). The simulated attosec- 



ond pulses do not have compact support, but exhibit a 
sharp Gaussian decay (See Fig[21 where the two original 
pulse polarizations are plotted in a logarithmic scale). 
Figure [3] shows the temporal profiles of the two polariza- 
tions of the simulated pulse and their reconstruction. As 
can be seen, their agreement is excellent. The errors are 
typically of the order of Err{T) ~ 0.05. 
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FIG. 3: Comparison between original (Amplitude - blue, 
Phase - black) and retrieved (Amplitude -dashed red, Phase - 
dash-dot green) pulses, for the two polarizations of the HHG 
pulse generated from the CO2 molecule at an alignment of 
58° (between the molecule symmetry axis and the IR polar- 
ization) and IR laser intensity of Q.Qlau at 800nm. The lower 
figure shows the radiation in the IR polarization, and the up- 
per shows the perpendicular polarization. The assumed pulse 
duration is 415 attoseconds. 

There are many experimental ways to implement the 
method developed here. So far we have discussed the 
spectrally dependent polarization components of HHG 
which can be measured, using an XUV polarizer, up to 
a sign ambiguity of the phase difference. The sign am- 
biguity can be resolved by means of a wave plate. For 
many cases this is not needed as by continuity the ambi- 
guity is the same for all uj, which means an overall time 



direction ambiguity. Another way to use the same idea 
is to generate two XUV pulses from different sources, ei- 
ther spatially [H, or by using mixed gases . By 
measuring the spectrum of each source alone and the 
interference between them, one can use our method to 
reconstruct the spectral phase. The spatial case can be 
generalized to spatio-tcrnporal measurements by lateral 
shearing interfcrometry [26| . 

To conclude, we have proposed and demonstrated a 
novel characterization method for attosecond pulses us- 
ing a vectorial phase retrieval algorithm. Our method 
presents a new class of solutions of phase retrieval prob- 
lems, applicable to many other fields, such as lensless 
imaging and optical spectroscopy. The main strength of 
the method lies in the fact that it removes the require- 
ment for a nonlinear interaction which is currently a ma- 
jor limiting factor in our ability to resolve many attosec- 
ond processes. Extending our approach to characterize 
more complex electron dynamics proposes a new scheme 
of time resolved measurements where attosecond-scale 
phenomena can be observed using linear, time-stationary, 
apparatus. 
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